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INTRODUCTION 


Proton- induced radioactivity of scintillation detectors is an 
important source of background- counting rate whenever these detectors 
are used in a space environment. In addition to the ever-present 
primary cosmic- rays, for spacecraft in Earth orbit, the trapped proton 
radiation adds considerably to the total proton flux. For orbits of low 
altitude and medium inclination, the largest fraction of the proton flux' 
is received in the region of the South Atlantic Anomaly. Figure 1 shows 
the flux contours for protons with energies greater -than, 5 MeV at an 
altitude of 600 kilometers (Ref. 1). Figure 2 gives the proton energy 
spectrum averaged over 24 hours for an orbit at an altitude of 720 
kilometers and an inclination of 3 0 degrees (Ref. 2). Figure 3 gives 
the primary cosmic- ray spectrum near the Earth. 

Radioactivity induced in spacecraft scintillation detectors was 
first recognized and analyzed by Peterson on his OSG I experiment 
(Ref. 3). The analysis of data from the Rus sian-gamma- ray experiments 
on the Cosmos-135 and Cosmos 163 spacecraft has included corrections 
for the ' radioactivity induced by passage through the region of the South 
Atlantic Anomaly, although the details of these corrections have not 
been given (Ref, 4), Dyer and Morfill (Ref. 5} have calculated the 
induced radioactivity in CsI(Ti) by the same general method used in 
.the present report and have compared these calculations to radioactivity 
induced by 155-MeV protons. This energy is representative of that of 
nuclear-interacting trapped protons. Their measurements are used in 
a study of the anticipated background radiation on the UK- 5 celestial 
gamma-ray experiment. The data contained within the present report 
has been used to determine the effect of cosmic-ray-induced radio- 
activity on diffuse celestial gamma- ray measurements (Ref. 6). 
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To deduce the amount of radioactivity and its energy and time 
decay spectrum induced in scintillation detectors by protons, both an 
analytical approach and an experimental approach, were undertaken. 
The analytical approach consisted of calculations of spallation yields 
and a study of the nuclear properties of the spallation products. The 
experimental approach was to irradiate a scintillation detector with 
high-energy protons and directly observe the counts resulting from the 
decay of induced radioactivity. Both approaches are limited in their 
estimation of the radioactivity induced in a true space environment. 
However, within errors, both estimates were in agreement and may 
prove useful in design and data reduction applications for space -borne 
experiments using scintillation detectors. 

Portions of this report were published in a paper entitled 
"Cosmic-Ray Effects on Diffuse Gamma- Ray Measurements" (Ref. 6) 
and were presented at the Joint Meeting of the. American Physical 
Society, Cosmic Physics Division and the American Astronomical 
Society, High Energy .Astrophysics Division in San Juan, Puerto Rico, 
on December 3 , 1971 . 
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CALCULATIONS 

/ 

Calculations of spallation yields were made for protons inter- 
acting with sodium iodide at energies above 100 MeV using the semi- 
empirical formulae derived by Rudstam (Ref. 7); These formulae give 
the cross section for spallation production within a factor of three over 
a wide range of atomic mass numbers and incident proton energies. 

The output of the program used to calculate the spallation production 

23 127 

cross sections for protons incident on Na and I targets is given in 
Appendix A; Each output group represents a different product element, 
and a matrix of cross sections is presented as a function of isotope 
number and interaction energy. Several qualitative features are 
apparent from the data. These features (listed, below) .are direct 
results of the basic Rudstam formulae (Ref. 7): 

• .The peak cross section for each element is close to 
to the valley of beta stability for that element. 

® At higher energies, the isotope production distribution 
is broader (cross sections fall off less rapidly) than 
at lower energies. 

e The element production falls off exponentially with 
' the difference in atomic number from the target 
element. 

® The above fall- off is much more rapid at Lower 
interaction energies. 

The output data are presented in pictorial form in Figures 4, 

5, and 6. Each of these figures represents a portion of the chart of 
nuclides (Ref. 8) in which the isotopes shown have a production cross 
section greater than that indicated on the figure at various energies. 
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’At an incident energy of 1 GeV, there are' 83 spallation nuclides 
produced in sodium iodide with yield cross sections greater than 
10 millibarns (1 millibarn = 1 0~ 27 cm 2 ) and 208 nuclides with yield cross 
sections greater than 1 millibarn. About 70 percent of these nuclides 
will either be beta- unstable or will have radioactive isomeric states 
with half-lives greater than 10 microseconds. (It should be noted that 
the total calculated spallation yield per unit crystal thickness for 
sodium iodide is'only 9 percent less than that of cesium iodide so that 
the results obtained would apply to either crystal. ) At proton energies 
other than 1 GeV, the cross section distribution of the spallation products 
changes greatly but the sum of the cross sections remains nearly 
constant, being approximately equal to the nuclear geometrical cross 
section of the target nucleus, 1, 260 millibarns for iodine and 410 millibari 
for sodium. 

-4* . 

A rough calculation may be made of the induced counting rate 
in a space craft- borne sodium iodide detector. The rate of detected 
counts per unit of isotropic detector area, F c , is given by 

F c = fi f r ft m c F p (1) 

isotropic flux of particles capable of inducing radioactivity 

-fraction .of F p which undergoes nuclear interactions 
within the detector 

fraction of spallation products which are radioactive 

fraction of. decays occurring within an allowable 
observation time 

average number of counts observed within the detector 
for each radioactive nucleus produced. 

Each term will now be considered in more detail. 


where 
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The term f| in Equation 1 may be calculated for thin targets 
composed of j different nuclear species from the formula 


% - 2 

j 


N o-j X j 


n 


where 



Avagadro's number 

cross section for a cosmic- ray spallation interaction 
with a target nucleus of atomic weight, A. 

J 


X- - mass thickness contribution of the target material, j, 

*) 

Assuming or. to be equal to the nuclear geometrical cross section and 
a typical detector thickness is = 25 gm-cm -2 , then f* = 0. 17. , This 

is sufficiently smaller than unity so that the thin target approximation 
is valid and indicates that cosmic- ray secondaries produced within the 
detector will have a minimal effect on the total spallation production. 

The terms f r , fj., and m c were derived from the nuclear decay 
schemes of over 100 spallation products, weighted according to their 
production cross section, as derived from the Rudstam formulae. It 
was found that the values were rather insensitive to the incident proton 
energy and that they are mainly determined by the nuclear energy level 
spacing of nuclei with mass numbers near that of the target nucleus. 

In determining ' fj., it was assumed that only those isotopes with half- 
lives in the range from 10 microseconds to 50 days would contribute 


to the observed counting rate. The values adopted are f r = 0.7, ft = 0.8, 
and m c = 1.6. The reason that m c is greater than unity is, of course, 
because daughter nuclei of many radioactive transitions are also 
radioactive. The relatively large value of ft results from the fact 





that the half- lives of most radioactive spallation products are in the 
minutes-to-days range. 

The half-life distribution of radioactive spallation products 
induced in an iodine target from 1 GeV protons is shown in Figure 7. 

In this figure, the sum of the cross sections of products within each 
logarithmic- spaced half- life interval is plotted as a function of the 
product half- life. 

Using the above values for the terms in Equation 1, the 
calculated counting rate within a 25 gm-cm" 2 detector (over all detected 
energies) relative to the inducing particle flux is 

F c 

~ — = 0.15 count per incident proton. 

* P 

This result includes only the effects of the residual target 

nuclei (spallation products) and not the smaller fragments which are 

ejected from the target nuclei. With the exception of Be 7 , which has 

a production cross section of approximately 10 millibarns at high energies 

all of the more abundant fragments that are produced are either very 

short-lived (n < 10 microseconds) or stable and thus do not contribute 
z ~ _ 

to the counting rate. 

The spectral distribution of the induced counting rate is not 
easily deduced from the nuclear decay schemes. However, it is noted 
that the spectrum of gamma rays from a large number of mixed 
fission products has an exponential form with an e- folding energy of 
0.9 MeV (Ref. 9). The average atomic weight of these products is close 
to that of iodine, and, thus, they have similar nuclear energy level 
spacings.' For the purpose of the present calculations, we will adopt a 
spectral distribution of the form 
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I 


dF 

c 

F p dE 


f p e * 


exp 



counts * MeV“ 1 


per incident. particle 


E e = 0. 9 MeV . 


The rather coarse assumption of a purely exponential energy- loss 
spectrum may introduce errors of up to a factor of two at any specific 
energy from 50 keV to 6 MeV, but the integrated counts in the spectrum 
are probably accurate to ±20 percent. Indeed, spectral features result- 
ing from electron capture decays of the more abundant radioactive 
products will certainly be present in the energy- loss spectrum. It 
should be noted that the observed energy loss will include the K shell 
binding energy of the parent isotope (33.2 keV for iodine). Listed in 
Table 1 are the expected spectral lines resulting from electron capture 
and internal transitions. 
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TABLE 1 . SPECTRAL LINES. FROM IODINE SPALLATION' PRODUCTS 


ISOTOPE 

DECAY* 

Xl / 2 


jl26 

EC 

13 d 

0.667 + 0.033 

j T 26 

EC 

13 d 

1.420 + 0.033 

r 

j 1 25 

EC 

60 d . 

0.035 +.0.033 

' i 124 

EC 

4.2 d 

0.603 + 0.033 

r 1 24 

EC 

. 4.2 d 

2.295 + 0.033 

i 124 

EC 

4.2 d 

1 .962 + 0.033 

jl 23 

EC 

13.3 h 

0.159 + 0.033 

-i 121 

EC 

2.1 h 

0.212 + 0.033 


18 IT 

58 d 

; » 

0.145 

Te 123m 

• IT 

117 d 

. 0.247 

Te 121nl 

-IT 

154 d 

0.293 

Te 121 ' 

EC 

' 17 d 

0.573 + 0.032 

Sb 122m 

IT 

4.2 m r 

0.026 

Sb 122 ™ 

IT 

530 usee 

0.075 

Sb 122m 

IT 

1 .8 usee 

0.061 


*EC = ELECTRON CAPTURE 
IT = INTERNAL TRANSITION 
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EXPERIMENT 


In an effort to directly measure the radioactivity induced in a 
scintillation detector, an accelerator experiment was performed at the 
NASA Space Radiation Effects Laboratory, Newport News, Virginia. 

A sodium iodide detector, 1.7 gm-cm" 2 thick, was irradiated with a ■ 
low flux of 600 MeV protons over a 10- second period. The integrated 
irradiation flux, 7 X 10 1 ® protons in a 5- centimeter- diameter sodium 
iodide detector, was measured using the known cross section for the 
I* 2 ^ (p, p 4n) X 123 reaction as a beam monitor (Ref. 10). The number 
of 1*23 nuc i e i produced were derived from the 0. 192- MeV line observed 
in the decay spectrum at a known time after the irradiation. Continuum 
and background radiation corrections were made to the data in this 
derivation; however, its accuracy is estimated at only ±30 percent 

i 

because of the magnitude of the corrections and' the uncertainty of the 
12 "? 

I cross section. 

Counting rate spectra were taken at various times after the 
irradiation. Several of these spectra are presented in Appendix B. 

A superposition of five of these spectra is shown in Figure 8. In 
general, the spectra show exponential continua with several discrete 
peaks superimposed. Most of the peaks from electron capture and 
internal transition decays of radioactive isotopes of iodine, tellurium, 
arid antimony in Table 1 could be identified in the spectra. The. 
continuous nature of the spectra arises from the large number of 
radioactive products involved and the fact that most transitions include 
a non- discrete beta transition. 

Intensity as a function of time after irradiation plots were made 
for various energy ranges from 0. 05 MeV to 3. 0 MeV and from 3.5 
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COUNTS (MeV-sec) 




hours to 50 days after irradiation. These data are presented in 
Appendix C, and a superposition of the decay curves is shown in 
Figure 9. Counts prior to 3.5 hours could not be measured accurately 
because of high counting rates. For each energy interval, the total 
counts were integrated from 3.5 hours to 50 days following irradiation 
and then multiplied by a factor 1 . 8 to account for the undetected counts 
from 10 microseconds to 3.5 hours after irradiation. The above factor 
was derived from the half-life distribution of spallation products formed 
by protons and sodium iodide at 600 MeV, according to the Rudstam 
formulae. This factor was assumed to be constant over the detected 
energies. The data were also normalized to a 25 gm- cm" detector 
thickness. 

Both the calculated cosmic- ray- induced counting rate and the 
600- MeV, proton- induced counting rate are plotted in Figure 10. The 
errors in the measured data, as indicated by the error bars on the 
lowest energy datum,, represent an uncertainty of a factor of 1.5. This 
uncertainty arises mainly from the determination of the irradiation flux 
and, thus, represents an absolute error applicable to each measured ' 
value. The relative error between the measured data is estimated to 
be ±20 percent. The error on the calculated data is estimated to be 
a factor of two, as indicated by the error bars in Figure 10. This 
factor includes both absolute and relative errors; the relative errors 

t 

arise from the uncertainty of the spectrum. 

The spectrum of the experimental counting rate is slightly softer 
than the assumed spectral shape; the data above 0.25 MeV can be best 
fitted by an exponential spectrum with E e 52 0.6 MeV. However, at 
each energy, the experimental and calculated induced counting rates 
agree to within the indicated errors. The energy independence of 
high-energy spallation radioactivity is supported by the measurements 
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of Barbier (R'ef. 11). Those measurements, indicate that uranium 
irradiated by 6 00 MeV protons produces radioactivity similar to that 
produced at an irradiation energy of 1’8 GeV. 
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SUMMARY -AND CONCLUSIONS 


The present study indicates that proton- induced radioactivity 
* \ 
will be a major source of background radiation for low-level, space- 

borne detectors in the hard X-ray and gamma- ray regions. Approxi- 
mately one gamma- ray will be produced after each nuclear interaction. 
The decay of the induced radioactivity will roughly follow a 1/t rate, 
and the energy spectrum is expected to be exponential with an e-folding 
energy of ~ 0. 9 MeV. Proton irradiation data at 600 MeV are in 
agreement with the above results. 

More precise calculations may be made using experimental 
cross section data and detailed decay schemes. The actual proton 
transport within heavy spacecraft will be needed for these calculations, 
along with the proton dose history. Such calculations are now in 
progress at the Space Sciences Laboratory /Marshall Space Flight Center. 
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APPENDIX A. Nal SPALLATION YIELDS 
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A- 1 1 
A- 12 
A- 13 
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A- 15 
A- 16 
A- 17 
A- 18 
A- 19 


SPALLATION YIELDS FOR Na ■ TARGET 
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APPENDIX B. 600-MeV PROTON IRRADIATION 
SPECTRA OF Naim) 

Pulse-height spectra were taken using standard commercial 
laboratory nuclear instrumentation. These spectra were recorded on 
paper tape and digital magnetic tape for later conversion into punched 
paper cards. The data on these cards were used in a computer program 
which normalized the data with respect to energy and accumulation time. 
The computer- generated histogram plots are reproduced in this appendix. 
The large peak at the end of several spectra are instrumental artifacts 
caused by amplifier saturation. Several discrete peaks in each spectrum 
are identified at the top of each spectrum. Refer to Table 1 for the 
energy and half-lives of these isotopes. 
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LOCATION* 

TIME AFTER 
IRRADIATION 
(Days) 

ACCUMULATION 
• TIME ' 
(sec) 

Background 

SREL 

Pre- Irradiation 

2, 000 

Background 

MSFC 

Pre- Irradiation 

40, 000 

1 

SREL 

0.29 

100 

2 

SREL 

0.41 

100 

3 

SREL 

0.60 

100 

4 

MSFC 

1.63 

1,000 

5 

MSFC 

2.5 

1, 000 

6 

MSFC 

5.4 

1, 000 

7 

MSFC 

7.0 

2, 000 

8 

MSFC 

15 

4, 000 

9 

MSF'C 

68 

80, 000 

10 

MSFC 

113 

82/000 


*SREL = Space Radiation Effects Laboratory, Hampton, Virginia 
MSFC = Space Sciences Laboratory, Marshall Space Flight Center 
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FIGURE B-l. BACKGROUND SPECTRUM, SREL MEASUREMENTS 
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FIGURE B-2. background spectrum, msfc measurements 
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FIGURE B-3 . SPECTRUM NUMBER ] 
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FIGURE B-4. SPECTRUM NUMBER 2 
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FIGURE B-5. SPECTRUM NUMBER 3 
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FIGURE B-6. SPECTRUM NUM8ER 4 
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FIGURE B-7. SPECTRUM NUMBER 5 
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FIGURE B-8. SPECTRUM NUMBER 6 
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FIGURE B- 10. .SPECTRUM NUMBER 8 
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FIGURE B-11. SPECTRUM NUMBER 9 


B-l 3 


COUNTS/ keV/sec 



E(keV) 

FIGURE B-12. SPECTRUM NUMBER 10 
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APPENDIX C. TIME DECAY CURVES FOR 600 MeV PROTON 

IRRADIATED NalfTA) 

The counting rate as a function of time after irradiation is 
shown for each of five energy ranges. Generally; the decay follows 
a 1/t function which flattens prior to the 0. 1 day and becomes steeper 

after 1 0 2 days. Counting rates in the 2- to 3- MeV energy range could 

<> 

not be taken .after 10 days after irradiation because of experimental 
diffi culties . 
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FIGURE C-l. COUNTING RATE AS A FUNCTION OF TIME AFTER 
IRRADIATION, RANGE 0.15 TO 0.25. MeV 
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FIGURE C-4; COUNTING RATE AS A FUNCTION OF TIME AFTER 
IRRADIATION, RANGE 0.6 TO 1.0 MeV 
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FIGURE C-5. COUNTING RATE AS A FI 
IRRADIATION, RANGE 1 
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FIGURE C-6. COUNTING RATE AS A FUNCTION OF TIME AFTER 
IRRADIATION, RANGE 2 TO 3 MeV 



